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The effect of the higher-energy 2nd resonance and the associated adiabatic-to-nonadiabatic transi- 
tion on neutrino propagation in solar matter is presented. For WIMP-annihilation neutrinos injected 
with energies in the "sweet region" between 300 MeV and 10 GeV at the Sun's center, a significant 
and revealing dependence on the neutrino mass hierarchy and the mixing angle #13 down to 0.5° 
is found in the flavor ratios arriving at Earth. In addition, the amplification of flavor ratios in the 
sweet region allows a better discrimination among possible annihilation modes of the solar dark 
matter. Under mild assumptions on WIMP properties, it is estimated that 200 neutrino events 
in the sweet region would be required for inferences of #13, the mass hierarchy, and the dominant 
WIMP annihilation mode. Future large-volume, low-energy neutrino detectors are likely needed if 
the measurement is to be made. 
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I. INTRODUCTION 

The annihilation of weakly interacting massive parti- 
cles (WIMPs), trapped in the Sun's core, is expected 
to produce fluxes of neutrinos and antineutrinos of all 
active flavors at energies well above the MeV scale of 
solar- fusion t/ e 's. These WIMP-annihilation neutrinos 
encounter the 2nd matter resonance at higher energy 
E% ~ |£m§ 2 | cos2<9i 3 /[2 V e (0)] - 0.2 GeV. We find that a 
small leptonic mixing angle 6*13 will inject a tell-tale sig- 
nal in the flavor spectrum at Earth in the energy region 
0.3-10 GeV, much as the 1st resonance at lower energy 
E l R ~ < 5TO^ 1 cos26»i 2 /[2V e (0)cos 2 6ii3] - 1.8 MeV is ex- 
pected to do in future measurements of the solar- fusion v e 
spectrum. This signal implies an experimental sensitivity 
to #13, presently constrained by CHOOZ data to be below 
12°, down to about half a degree. In the above expres- 
sions, 5m 2 jk = m 2 — m|, and V e (0) ~7x 10~ 12 eV is the 
matter potential resulting from the electron density at 
the Sun's core; the superscripts h and I denote the higher- 
energy 2nd resonance and the lower-energy 1st resonance, 
respectively. The signal comes from the adiabatic-to- 
nonadiabatic transition at the 2nd resonance. This effect 
is the higher-energy boundary of the "bathtub" spectral 
shape, well-known and well-described for the 1st solar 
resonance, and elucidated for the 2nd solar resonance in 
Figs. 3 and 4 of Ref. [l| as well as herein. 

The Borexino experiment Q in progress and the SNO+ 
experiment [3j in the construction stage are likely to 
measure the 7 Be and pep monochromatic solar neutri- 
nos at energies of 0.861 MeV and 1.442 MeV, respec- 
tively. These measurements should reveal the lower- 
energy boundary of the bathtub profile from the 1st res- 
onance. The adiabatic-to-nonadiabatic transition at the 
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2nd resonance offers a potentially much more striking sig- 
nal than the lower-energy resonance, as we will show in 
the present work. The heuristic reason is that MSW reso- 
nant enhancement of the small 6*13 to the MSW resonant 
value of 45° is a much larger effect than the enhance- 
ment of #12 ~ 32° to the resonant 45° value. In fact, as 
we will demonstrate, the adiabatic-to-nonadiabatic tran- 
sition at the 2nd resonance is extremely sensitive to the 
value of small #13: the resonance is completely nonadi- 
abatic and ignorable for #13 = 0, whereas the adiabatic 
component is significantly amplified for nonzero #13, even 
as small as 0.5°. In addition, the 2nd resonance occurs 
in the neutrino sector if the neutrino mass hierarchy is 
normal (i.e., 5m^ 2 > 0), but in the antineutrino sector if 
the mass hierarchy is inverted (i.e., Sm^ 2 < 0). WIMP 
annihilation in the Sun is expected to produce neutrinos 
and antineutrinos in equal numbers. The different ab- 
sorption and scattering cross-sections of neutrinos versus 
antineutrinos in terrestrial detectors then allows possible 
discrimination between the two. It thus becomes appar- 
ent that the 2nd resonance provides the potential to re- 
veal not only the value of #13, but also the neutrino mass 
hierarchy codified as sign(<5m 2 2 ). 

As in the case of solar-fusion neutrinos, the total 
neutrino spectrum remains unaffected by the resonance 
structure [4J . It is rather the distribution of the neutri- 
nos over the three flavors that depends on the resonance 
physics. For this reason, we focus our attention on the 
individual flavor spectra (three flavors each for neutri- 
nos and antineutrinos) that experiments at Earth can, in 
principle, detect. 

While much theoretical effort has gone into studies of 
"indirect detection" of solar WIMPs by identification of 
their high-energy neutrino flux in large terrestrial de- 
tectors [a], little has been done [3, [(| to elucidate the 
possibilities at lower energies, E < 10 GeV. The rea- 
sons are clear: most of the neutrino flux from WIMP 
annihilation is expected to populate the higher-energy re- 
gion, the atmospheric background falls as ~ E~ 3 above 
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10 GeV Q, and the detection cross section for a neutrino 
grows as ~ E above GeV neutrino energies. However, 
the substantial discovery potential that hides in the solar 
WIMP-annihilation neutrino data below 10 GeV compels 
us to promote the associated lower-energy physics. In the 
next section, we present an overview of this rich neutrino 
physics at lower energies. Details, drawn mainly from 
Ref. [l|, are presented in subsequent sections. 



II. PRELIMINARIES AND OVERVIEW 



A statistical average over the oscillation phase <f> = 
dm^ k L/2E effectively results from the uncertainties in 
the baseline L and energy E of solar neutrinos. The pro- 
duction flavor ratios w a (E) in the Sun's core and the 
terrestrial flavor ratios W a (E) for neutrinos are then re- 
lated by 0: 



W» ) =UPUl(r = 0) 




(1) 



Here, (U m ) a j = \(U m ) a j\ 2 is a classical probability ma- 
trix, constructed from the mixing matrix in matter U m . 
The transposed matrix (r) transforms the production 
flavor fluxes to the propagating mass-state fluxes. The 
argument r — in Eq. ([1]) reminds us that WIMP anni- 
hilation occurs in the solar core, and so it is the matter 
density at the Sun's center that determines U^(0). 

Possible nonadiabatic transitions between the 
effective-mass states at the 1st and 2nd resonances and 
other solar-matter effects are described by the level- 
crossing probability matrix P [8l4Icj] . Finally, the matrix 
determined from the usual vacuum mixing 
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matrix U, transforms the mass states back to detectable 
flavor fluxes at Earth. Explicit expressions for U m and P 
in terms of neutrino energy E, solar-model coefficients, 
and neutrino parameters, as well as the analogous result 
for antineutrinos can be found in Ref. along with 
many other relevant details. Some related references 
that study flavor issues for solar neutrinos arising from 
WIMP annihilation are listed in Ref. fill ]. 

Inspection of Eq. (fT]) reveals that structures in the 
flavor spectra can potentially arise from three energy- 
dependent sources: the initial flavor ratios w a (E), the 
mixing matrix in matter U m (r = 0,E), and the jump- 
probability matrix P(E). The vacuum mixing matrix U , 
and thus U, are each independent of the neutrino en- 
ergy, and so U does not contribute to energy-dependent 
features in the terrestrial flavor spectra. 

The first of the above sources for features in the terres- 
trial flavor spectra, namely w a , is WIMP-model depen- 
dent. A description of this energy dependence requires 
knowledge of the WIMP— > v a annihilation chains. Three 
possible decay chains are commonly invoked. They are 
neutrino production via decay of the intermediate states 
W + W~, bb, and t + t~. The latter two chains lead to 
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FIG. 1: Shown are the relative neutrino flavor spectra re- 
sulting from annihilation of WIMPs with 100 GeV mass into 
predominantly t + t~. This plot is based on the numerical 
results obtained by Cirelli et al. [12| ■ The shaded area repre- 
sents the "sweet region" defined in the text. Note the nearly 
constant values of the relative fluxes inside this region. 



softer neutrino spectra than does the former. Calcula- 
tions 12j reveal that for the bb and t + t~ annihilation 
chains, the neutrino flavor ratios at production are slowly 
varying in the low-energy region |l3[. For the W + W~ 
chain, these flavor ratios are more energy-dependent [3] . 
We show the nearly energy-independent neutrino flavor 
spectra for Mwimp = 100 GeV annihilation via the t + t~ 
mode in Fig. [TJ For a 100 GeV WIMP mass, the region 
of nearly constant flavor ratios lies below 20-30 GeV [15] ■ 
For specificity in the rest of this paper, we will continue 
to focus on the bb and (especially) t + t~ modes, as well 
as on a WIMP mass of order 100 GeV. If the dominant 
annihilation mode of the WIMP is W + W~ rather than 
t + t~ and/or bb, then we expect that the resonant flavor 



change that we describe below can still be extracted from 
data, but further efforts would be required to isolate the 
resonant features from the non-resonant energy-varying 
flavor ratios. 

The second potential source of energy dependence, the 
matrix U m , can also be tamed: for neutrino energies 
sufficiently above the 2nd resonance, U m approaches an 
E- independent constant matrix (see e.g., Ref. [l]). For 
this reason, we will consider only energies well above 
E^ ~ 0.3 GeV, so that U m is effectively constant. 

Some energy dependence will also be introduced in the 
flavor ratios by the adiabatic-to-nonadiabatic transition 
of the 1st resonance. This occurs at energies at and above 
E l NA ~ 10 GeV PJ. As can be seen in Figs. 3 and 4 
of Ref. [l|, the effects of the 1st resonance are mild; to 
ignore these effects would introduce little uncertainty into 
our calculation, and the effects could be included in a 
more complete analysis. However, we will simply restrict 
our study to energies below E l NA to shield our analysis 
from "contamination" due to the 1st resonance. Thus, 
we arrive at, and define, a "sweet region" in energy for 
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our analysis: 

0.3 GeV < E < 10 GeV "SWEET REGION." (2) 

The results presented below will hold within this special 
region of energy. 

The remaining and most important source of struc- 
ture in the flavor spectra arises from the neutrino level- 
crossing probabilities at the resonances, which are de- 
scribed by the jump matrix P. The transition from the 
adiabatic (i.e., no level-crossing) to the nonadiabatic (i.e., 
complete level-crossing) regime at high energies can leave 
a dramatic imprint in the terrestrial flavor-flux ratios. 

In Ref. [l|, the energy E^ A characterizing the onset 
of this adiabatic-to-nonadiabatic transition was defined 
implicitly by setting P^(Ej^ A ) in the crossing-probability 
matrix to e~ 3 ~ 5%. Numerically, this onset energy is 
E^ A ~ 750 sin 2 #i 3 GeV. (An expression for the crossing 
probability is given below in Eq. (fl3|) .) The 2nd reso- 
nance occurs in the neutrino sector for the normal hier- 
archy with Sm^2 > 0, and in the antineutrino sector for 
the inverted hierarchy with <5m 2 2 < 0. 

The above discussion shows that we may identify any 
observable energy dependence of the flavor spectra in 
the sweet region as due to the adiabatic-to-nonadiabatic 
transition at E^ A . In turn, the energy dependence in 
this sweet region will implicate the value of a small but 
nonzero #13 as well as the neutrino mass hierarchy. In 
particular, at small 0x3 the nonadiabatic jump probabil- 
ity is exp(— r# 2 3 ). The "adiabaticity parameter" T can 
be written in terms of the length scale Xq = \4- In A^l -1 
for the change in the electron density iV e at the reso- 
nance region (corresponding to one e-folding in density 
for an exponential profile like that of the Sun) , and the 
neutrino oscillation length in vacuum A„ = An E/Sm 2 as 
r = 47r 2 Xq/X v . The important feature for our purposes 
is that r scales as 1/E and is therefore large, roughly 
300-10 000, across the low energies of our sweet region. 
It is this large size of T that produces an observable sig- 
nal even for extremely small #13. Qualitatively, we expect 
sensitivity to #13 in the sweet region for #13 as small as 
VVT > 10~ 2 - 0.5°. In Sec. [TV] we show quantitatively 
that this is indeed the case. 

There is also a terrestrial source of energy dependence 
in the flavor ratios. This arises from matter effects for 
neutrinos and antineutrinos traversing the Earth. The 
effect has been well worked out in publications [l6| ; it is 
somewhat complicated (best addressed with a numerical 
code), and we will not include it in this paper. We do 
give here some general remarks jlj about the Earth mat- 
ter effect that are relevant for our present purposes. The 
resonant energies in Earth related to the solar scale 5m|i 
are below the sweet region (at ~ 25 MeV and ~ 100 MeV 
for the Earth's core and mantle, respectively), while the 
resonant energies in Earth related to the atmospheric- 
scale 5m| 2 fall right in the sweet region, at ~ 2 GeV and 
~ 10 GeV for the Earth's core and mantle, respectively. 
As with the solar resonance, for a normal mass hierarchy 
the relevant resonance in Earth occurs in the neutrino but 



not the antineutrino sector. And again as with the so- 
lar resonance, for an inverted mass hierarchy the relevant 
resonance in Earth occurs in the antineutrino but not the 
neutrino sector. The resonant amplification can be quite 
large, but the effect can be mitigated by a resonant os- 
cillation wavelength scaling like 1/613. The requirement 
that at least a quarter of the oscillation length must lie 
within the Earth to "feel" the matter effects leads to a 
nadir-dependent condition on #13. For a neutrino with 
zero nadir angle, 813 > 0.5° (E /GeV) is required to feel 
the matter; for other nadir angles, larger values of 6*i 3 are 
required to feel the matter. In a more complete analysis 
involving real data, terrestrial-matter effects should then 
no longer be ignored. 

It is worth remarking here what would change if the 
WIMP mass were, say, a TeV rather than the 100 GeV 
we have assumed [13|]. Surely, the available phase space 
for the produced neutrinos is increased. Moreover, in 
Ref. [l2j it is shown that the upper end of the neutrino 
energy spectrum scales linearly with the WIMP mass, 
tending to enlarge the energy region in which the flavor 
ratios are constant. However, the sweet region remains 
between 0.3 and 10 GeV because it is determined solely 
by the effects of solar matter on neutrino propagation. 
It follows that the phase space of the sweet region rela- 
tive to the enlarged total phase space now represents a 
smaller fraction, so that less sweet-region neutrinos are 
available. More importantly, according to Eq. © the ex- 
pected neutrino flux from WIMP annihilation in the Sun 
scales as M^ 2 MP , and so would be down by an additional 
factor of 100 for annihilation of TeV WIMPs compared 
to 100 GeV WIMPs. 

One must inevitably ask what event rate at Earth 
might be expected from solar WIMP annihilation to 
(anti) neutrinos in the sweet region, and what rate is 
needed to observe structure in the terrestrial flavor spec- 
tra. We attend to first of these questions in the next 
section, and discuss the requirement for statistical signif- 
icance in our final Sec. IVII 



III. (ANTI) NEUTRINO FLUX FROM SOLAR 
WIMP ANNIHILATION 

Theory suggests that the age of the Sun exceeds the 
equilibration time between solar capture of WIMPs and 
their subsequent annihilation in the Sun (l7j . Conse- 
quently, the WIMP annihilation rate is given by half of 
the WIMP capture rate, where the "half" just reflects 
the fact that it takes two captures to enable one two- 
body annihilation. The WIMP capture rate by the Sun 
is given in Ref. fl% : 

C q ~ 1.0 x 10 25 .Fs" 1 , (3) 
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where 
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/ 100 GeV \ 2 

V -^WIMP / 



(4) 

is a fiducial factor. Here, ctwimp and Mwimp are the un- 
known WIMP-nucleon scattering cross section and mass, 
while pwimp and t»wiMP denote the density and rms ve- 
locity of the local WIMP population. Each parenthetical 
fraction displays typical values, except that the fiducial 
value shown for ctwimp is the present upper limit for the 
spin-dependent cross section [Tj| 0j| ; it must therefore be 
viewed as optimistic. The inverse quadratic dependence 
of the rate on the WIMP mass is easily understood: one 
factor arise from the conversion of WIMP mass density 
Pwimp to number density, and the other factor is kine- 
matic, reducing the capture efficiency when the beam 
and target masses are mismatched. The fiducial values 
for today's local WIMP density and rms velocity serve as 
numerical guidelines, but in fact the integrated history 
of WIMP capture by the Sun over large look-back times 
may yield values that do not adhere to these guidelines. 

The WIMP-nucleon cross section in Eq. (Q} above is 
the sum of a spin- independent cross section and a spin- 
dependent cross section, each averaged over the target 
matter in the Sun's core. Experiments on Earth that 
search directly for WIMPs typically use target materials 
with large nucleon number, and so are better suited to 
limit (or detect) the spin-independent cross section (20j . 
Direct-search limits on the spin-dependent cross section 
are five orders of magnitude weaker than the limits on 
the spin-independent cross section. Since the solar mate- 
rial is mainly hydrogen, solar capture of WIMPs is very 
sensitive to the spin-dependent cross section. 

An experimental measurement of the solar neutrino 
flux from WIMP annihilation would bypass the theoret- 
ical uncertainties just listed. So far, experiments have 
yielded only upper limits on this flux inferred from final- 
state muons with energies above a GeV [18|, ll9j . It is in- 
teresting to note that the bounds on the spin-dependent 
WIMP cross section inferred from these experimental 
neutrino-flux constraints are stronger than the bounds 
coming from direct searches for WIMPs plj . 

We continue by considering the mean multiplicity of 
neutrinos, £„, produced per annihilating WIMP (not 
WIMP pair), a quantity we anticipate to be of order 
unity. We can then estimate the total neutrino flux at 
Earth 
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The background from "atmospheric neutrinos," i.e., the 
neutrinos resulting from the decay of charged pions pro- 
duced by cosmic-ray interactions in our atmosphere, 
has been calculated by many groups with convergent 
results. A typical energy spectrum in the GeV re- 
gion can be found in Ref. Q. It approximately obeys 



1.2 x 10 6 (Jv) 275 (GeV • yr • cm 2 )" 1 for + the 
all-flavor flux would be about three times larger. The 
integrated all-flavor flux above a GeV is then about 
2 x 10 6 /yr cm 2 , with the spectral-index factor nearly 
compensating the flavor factor. This background flux is 
roughly twenty times the maximally allowed solar WIMP 
neutrino flux. Moreover, the background flux has a fla- 
vor content that varies with direction on the sky. The 
down-coming atmospheric neutrinos will show the 1 : 2 
v e to 1/^ flavor ratio characteristic of the complete pion 
decay chain, as the muon decay length at lower energies 
is shorter than its atmospheric height, and the neutrino 
pathlength is shorter than its vacuum oscillation length 
L osc = 4nE/Sm 2 ~ 1000 (E/GeV) km. In contrast, the 
upcoming neutrinos will show a 1:1 flavor ratio, since half 
of the v^s will have oscillated into iVs. On the bright 
side, the solar fraction of solid angle on the sky is quite 
small, approximately [n (R & / AU) 2 ] / An ~ 5.4 x 10~ 6 . So 
one may hope that a cut favoring the direction to the 
Sun would greatly increase the signal-to-background ra- 
tio. However, the mean scattering angle for neutrinos 
below 10 GeV is large, (9 scatt ~ 20° y^TOGcV/E. The 24- 
hour angular modulation expected in the solar signal will 
help reduce the unmodulated atmospheric background. 
In addition, there has been some recent discussion of 
possible very large detectors [23| , and some development 
in directional reconstruction of neutrinos at lower ener- 
gies [H], [H[ . There have also been recent studies of using 
event topologies to achieve partial "statistical" separa- 
tion of neutrino and antineutrino data samples (26j . 

Multiplying the solar WIMP neutrino flux by (i) the 
neutrino-nucleon cross section, by (ii) the target number 
of nucleons given by 



N N = 6.0 x 10 
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megaton J 
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and by (iii) the fraction of incident neutrino flux in the 
sweet region A£(SR) from 0.3 to 10 GeV 



r Win 7? dN " 

_ JAE(SR) amrj d\nE 
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(7) 



one obtains the event rate within the sweet energy region. 
For a 100 GeV WIMP mass, inspection of the theoreti- 
cal neutrino spectra in Ref. [TJ[ suggests a value ~ 20% 
for /„, the fraction of neutrino flux in the sweet region 
between 0.3 and 10 GeV. As a fiducial event rate, we 
therefore take 



R = 140 



<?vn{E) 



target 



10~ 38 cm 2 / V megaton / \20% 



( u \ iuT 



yr 



(8) 



For some perspective on the 140 events per year, we may 
ask how many events per year are to be expected in a 
proton-decay experiment at a megaton detector. The 
present limits on the lifetime of protons to decay to var- 
ious modes are typically 10 35 yr. The nucleon number in 
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a megaton is 10 12 Na — 6 x 10 35 . So the expected event 
rate is 

ZlO^yA / Mdetector \ _ X , Q , 

-Rp+docay = 6 — yr . (9 

y 3 \ t p J \ megaton / 

We see that the rates for detection of solar neutrinos in 
the sweet region from WIMP annihilation and for de- 
tection for proton decay are comparable (although the 
backgrounds are different). 

IV. THEORETICAL FRAMEWORK 

To proceed, we need explicit expressions for the neu- 
trino flavor spectrum produced by WIMP annihilation 
in the Sun's core in the range 0.3-10 GeV. We begin 
by parameterizing the relative production fluxes w a . In 
many models, the dominant WIMP-annihilation chan- 
nels are bb, t + t~, and W + W~. The neutrino spec- 
trum from each channel has been calculated (see. e.g., 
Ref. [Hj]), but the branching ratios to these channels de- 
pend on the specific model. Fortunately, WIMP d ecay 
obeys w e (E) = w fl (E) to a very good approximation [12j . 
Since the normalization w e ± ± w T = 1 holds at any 
energy, the parametrization of the relative flavor spectra 
at production requires just one WIMP-model-dependent 
function /f{E). 

We define fp operationally via w e = = ^ — fp and 
w T = ~ + 2 fp. Positivity of the w a then implies the 
bound — | < fp < +|. As defined, fp is the deficit of 
w e or Wp from | at injection in the solar core. Since 
maximal mixing effects symmetrization of v„ and v r , it 
is also useful to view fp as the excess of uy = ± w T 
from |. To order fp, the ratios (w e /w T ) and (w^/wr) at 
injection are 1 — 9 fp + 36 f F . 

For comparative purposes later, it is also useful inter- 
pret fp in terms of the evolved flavor probabilities W a in 
the absence of matter effects. The flavor density matrix 
at injection, in the flavor basis, is given by 

PF = \t + f F [2\v T ){ VT \-\^){ Vtl \-\v e ){v e \] . (10) 

Using the tribimaximal mixing values to write this in the 
mass basis, and invoking phase averaging to remove the 
off-diagonal elements, one is left with just 

PF = \i + \f F • (ii) 

Then, the relative flavor probabilities W a are just 
( v a | pp | i> a ) . The results for vacuum transitions are 
W e = W e = | - § f F , and W„ = W T = W M = W T = 
g + \ fF- The vacuum value for the ratios W^/We and 
W^/We to order f 2 F are 1 + § f F + /|. 

Although actual WIMP properties are unknown, we 
may consider two of the aforementioned annihilation 
channels, bb and t + t~, as generic examples. We have 



TABLE I: Values of f F = | - w e = \ - to p = \{w T - |) in 
the sweet region for popular WIMP annihilation modes. 



WIMP annihilation mode fp (sweet region) 

t + t~ +0.182 ±0.001 

bb -0.09 ±0.01 

Higgs-like -0.07 ± 0.01 

W + W~ -0.029 ± 0.014 (E/GeV) 



argued that in our sweet region, the function fp (E) is re- 
duced to an energy-independent number fp for the t + t~ 
and bb annihilation modes of the solar WIMPs. Our fit- 
ted values for fp in these channels are —0.09 ± 0.01 and 
±0.182 ± 0.001, respectively. 

It is worth mentioning that in some WIMP models 
the annihilation to fermions proceeds through Higgs- 
like couplings. We then expect the bb mode to dom- 
inate, but also a branching fraction to t + t~ given by 
(m T /mi,) 2 /3 ~ 5%, where the | reflects the bb mode's 
color factor. In the sweet region, the /f's for the bb 
and t + t~ modes have opposite signs, so the average fF 
for the fermionic channel is smaller in magnitude than 
/f (bb) and /f(t + t~). However, the bb mode, and thus 
its concomitant fp value, dominate: an estimate employ- 
ing the injection spectra given in Ref. [l2| indeed yields 
/^(Higgs) ~ -0.07±0.01 in proximity to the -0.09 value 
for b b mode. We remark that the neutrino spectrum 
from the t + t~ mode is harder than that from the bb 
mode, and at energies above the sweet region, of no rel- 
evance to the present paper, the r + r~ mode becomes 
increasingly important and eventually comes to domi- 
nate. For the W + W~~ mode, fp behaves quite differ- 
ently: it rises nearly linearly with energy in the sweet 
region: f F (W+W-) = -0.029 ± 0.014 (E/GeV). The 
values for the various /f's presented in this paragraph 
have been estimated using the results in Ref. [12|; they 
are also collected in Table |U 

We note that for equal initial flavors (/i?=0), the re- 
sulting terrestrial flavor ratios are also democratic [l|, 
i.e., W e = — W T — 1/3. This is most easily seen by 
noting that the commutator in the density-matrix evolu- 
tion equation i -£ = [H, p] vanishes for p proportional to 
the identity matrix regardless of whether H is the vac- 
uum or matter Hamiltonian. Thus, it is only the order 
7-18% /p-valuc differences in the initial flavor spectrum 
that evolve nontrivially. If flavor differences were not am- 
plified by intervening matter effects, an event sample of 
about (n/15%) 2 would be necessary to yield an n-sigma 
statistical inference of this difference. However, we will 
see shortly that matter effects, quadratically sensitive to 
#i3 and to the neutrino mass hierarchy, may amplify the 
flavor difference considerably. 

For further progress, we also need explicit expressions 
for U, P, and U m appearing in Eq. (flj. In the sweet 
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TABLE II: Leading-order contributions to the components of the A and B flavor vectors in Eq. (|12[) . 

type e flavor flavor r flavor 

A™ 2 + 4(56*23 -1-85023 -l + 450 23 

-2-4^50 i2 + 4 50 23 1 + \/8<Wia + 856*23 + \/80i 3 cos 5 1 + y/E 80 12 - 12 69 23 - \/8 13 cos 6 

A™ -2 + 4^2 56>i2 +4 56/23 1 - ^50i2 - ^8 0i3 cos 5 1-^8 5012 + 78 013 COS 5 — 4 5023 



A™ -4V2 50 12 + 8 5023 \/8 50 12 + s/% 013 cos 5 \/850i2 — -\/80i3 cos5 — 85023 

B™ 

A 1 * 2 - 4 5023 -1 - 4 5023 -1 + 8 5023 

~B™ -4 + 4^2 5012 + 8 5023 2 - V8 56 12 + 4 50 23 - \/89i 3 cos 5 2 - \/8 50i2 - 12 50 23 + a/8 013 cos 5 



region, the level-crossing probability at the lower reso- 
nance is zero, as explained above. The fact that Eq. ([T]) 
is linear in the level-crossing probability matrix P allows 
one to write a simple equation for the flavor evolution, 



W a = \ + \[A a 



B a P?(E)]f F , 



(12) 



valid in the sweet region. The pre-factor of ^ in front 
of the square brackets is chosen for convenience. The 
flavor-indexed quantities A a and B a , a = e, fx, r, are de- 
termined by the neutrino mixing parameters at the Sun's 
core (where v e ~ z/3) and at Earth (vacuum values) [27| . 
Since the mixing parameters are energy independent in 
the sweet region, so too are A a and B a . The flavor- 
vectors B a and A a do depend on the neutrino mass hier- 
archy and on incident neutrino versus antineutrino. Im- 
portantly, the entire energy dependence in the evolution 
equation (|12j) is contained in the level-crossing probabil- 
ity Pc{E) at the 2nd resonance, given by the expression 



P?{E) = Q(E-E 



'R 



exp(-rsin 2 i3 ) - exp(-T) 



1 - cxp(-r) 



exp(-rsin 2 13 ) 



(13) 



In this expression, denotes the unit-step function, and 
0(Sm2i/\Sm^ 2 \) ~ 0.03 terms have been dropped. As 
mentioned in the overview section, the adiabaticity pa- 
rameter r may be written as 



T = 4tt 



2 ^0 



where 



A m = 



— lnJVe(r) 
ar 



(14) 



(15) 



is the scale of density change (the distance for an e- 
folding change in the solar density, equal to 6.6 x 10 4 km) 
and A„ = 47r E/Sm 32 is the oscillation length in vacuum. 
Consequently, T scales as 1/E, and we define a constant 
r via 



r 



r 

E 



where 



2.5 x 



\8ml 2 \ 



2.5 x 10- 3 eV" 



TcV. 



(17) 



for the higher-energy resonance in solar matter. With 
Eqs. ([li | -(frf |) at hand, it is apparent that T > 300. 
This large T value, together with the experimental input 
sin 2 #13 <C 1 and the fact that E > E^ in the sweet re- 
gion, has been used to obtain the second line of Eq. (fKl]) . 

Note that the only neutrino mixing parameter con- 
tained in Pc(E) is 6*i3, and it enters essentially squared in 
an exponent. It follows then, that the dominant energy 
dependence of the terrestrial flavor ratios in the sweet 
region is governed by #13, as advertised in the introduc- 
tion. Moreover, the sensitivity of flavor evolution through 
the adiabatic-to-nonadiabatic transition is exponentially 
sensitive to 6* 2 3 . 

From Eqs. (|13p and (fTB|) . we infer the onset energy for 
E^ A for the adiabatic-to-nonadiabatic transition to be 



E 



h _ 
NA = 



i f sin 2 e 13 



800 sin 2 9 



13 



|5m| 2 | 



2.5 x 10" 3 eV 2 



GeV. (18) 



Neutrinos with energy E^ < E < Ej^ A will experience 
adiabatic level-repulsion at the 2nd resonance, whereas 
neutrinos with E 3> E^ A will experience nonadiabatic 
level-crossing at the 2nd resonance. We remark that 
although it is true in the Sun that the adiabatic-to- 
nonadiabatic transition energy E^ A exceeds the resonant 
energy Ej>, this ordering is not in in general guaranteed. 
The ratio of the two energies is 



rph 
^NA 



2?r sin 2 6 



13 



3 cos 20i3 



V2G F N e 



V2G F X Q N e (0) 



(19) 

Numerically, the term in brackets is 1800 for the Sun. It 
is this fortuitously large value for the Sun that allows the 
adiabatic-to-nonadiabatic transition to probe $13 values 



( 16 ) all the way down to 9 13 ~ 1800)^ - l c 
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We proceed by defining the energy Ei/ 2 at which the 
crossing probability is one half. From Eqs. (TT5]) and p^|) 
we readily find that 

fsin 2 fl 13 „ R . 2 , ( \$mh\ \^, T 

Ei 11 — — ; — 7. — = 3.6 sin #13 tt lev . 

1 In 2 V2.5xfO- 3 eV 2 / 

(20) 

Let us also define the width of the change in crossing 
probability as AE = E + — E_ , where P£{E- ) = 1/e and 
P*{E+) = 1 - f/e. This yields 

A£ = (e-l)(ln2)£ 1/2 ~ 1.19 £ 1/2 . (21) 

It is seen that E 1 / 2 as well as AE are approximately 
equal to 4 sin 2 9 13 sa 4 6> 2 3 in units of TeV. Both E 1 / 2 
and AE are very sensitive probes of the unknown #13. 
The experimental upper bound sin #13 < 0.22 at 3<r es- 
tablishes that 813 ~ sin 6*i3 to better than 1% accuracy. 
Accordingly, we may invert Eq. (|2"0|) to obtain 

For #i3 in the phenomenologically interesting range 
0.6° < 6<i3 < 3°, Eq. (gDJ tells us that £i /2 lies within the 
sweet region. Then, the inverse Eq. (|2"2"j) establishes that 
a measurement of the adiabatic-to-nonadiabatic transi- 
tion energy, via changing flavor ratios, probes the #13 
mixing angle in this very interesting range from 0.6° to 
3° . Equation (|2"2"|) is the most important result of this pa- 
per. We show below that the flavor signal is also sensitive 
to the neutrino mass hierarchy. 

The ratio AE/E1/2 — 1.19 is independent of 6*13. Con- 
sequently, the number of events needed to establish i?i/ 2 , 
and thereby $13, is independent of 6*13 or E 1 / 2 . Of course, 
the necessary event number will depend on the magni- 
tude of flavor change over the transition region. Since 
AE ~ Ey 2 , we have that AE/E = AlnE ~ 1. From 
this we infer that the relevant event sample must span 
about a natural log in energy about Ei/ 2 to map out the 
transition region in detail, again independent of 8i 3 or 
Ei/2- However, to measure the gross change in flavor ra- 
tios across the transition, an event sample may be drawn 
from an energy region as large as the entire sweet region, 
0.3 to 10 GeV. The availability of this wide energy range 
for the data sample is fortunate, for at the relatively low 
energies of the sweet region, large uncertainties are ex- 
pected in the reconstruction of the neutrino energy from 
the measured charged lepton. 

The change in the flavor ratios across the adiabatic-to- 
nonadiabatic transition region of width AE centered on 
E1/2 is the signal. The magnitude of the flavor change 
across the transition region is determined by Eq. ([12"]) 
and depends on the parameters A a and B a . These pa- 
rameters in turn depend on the neutrino mixing angles, 
the mass hierarchy, and the particle type (neutrinos vs. 
antineutrinos). The two hierarchies and the neutrino ver- 
sus antineutrino type lead to four possibilities. To distin- 
guish among these possibilities, we employ the notation 



4™, A™, etc., where the bar denotes the antineutrino 
case and the superscripts NH and IH refer to the nor- 
mal and inverted hierarchy, respectively. For our present 
purposes, it is sufficient to employ the "tribimaximal" 
values |28| for the two large mixing angles, i.e., to set 
sin 2 8j 3 BM = 1/2 and sin 2 8 J 2 BM = 1/3. The phenomeno- 
logical accuracy of this approximation is very good, as is 
summarized in Ref. [29(. 

Although we employ the tribimaximal values for 832 
and 8\2 in what follows, we nevertheless exhibit in Ta- 
ble mi for possible future use, the first-order corrections 
to the A a and B a that result if the neutrino mixing an- 
gles deviate from their tribimaximal values. To this end, 
we have denoted the deviations from the tribimaximal 
case by 58 jk = 8jj~ — 8j BM . The expressions in Table UU 
result from a calculation paralleling that in Appendix D 
of Ref. [l|. The 3<r experimental uncertainties in the an- 
gles are [3(| 

- 10.4° < 58 32 < 8.6° , 
-4.8° < 6821 < 4.0° , 

0° < 0i3 < 12.9° . (23) 

Consequently, future data may mandate deviations from 
the tribimaximal-based A a and B a used in this work. 



V. RESULTS 

In the first subsection, we present results valid to all or- 
ders in fp. In particular, we provide figures for the t + t~ 
annihilation mode generated with the complete energy 
dependence of /f (E) over an energy range that includes 
the sweet region. In the second subsection, we present 
analytic formulas for the sweet region valid to order f F 
and compare (favorably) their validity with the more ex- 
act figures. 



A. Exact results 

The explicit expressions for the neutrino and antineu- 
trino flavor ratios with the normal hierarchy are 





1 


^(2-2P c h )f F 




</ H . = 


1 

3 


h^(-H-[l±V8. 


013 cos 5} P*)f F 


W NH = 

* v e 


1 

3" 


^(-2)/f 




w NH 

W »/T = 


1 
3 


~j(1tV88 13 cos 


S)f F - 
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The analogous expressions for the flavor ratios with the 
inverted hierarchy are 

= ~+0-/ F (25) 

Wl% = ±±±(V89 13 co S 5)f F 

Wf = I + i(2-4P c ")/ F 

W™ /T = ^ + ^(-l + [2 T V8e 13 cosS}P^f F . 

In these equations, the upper and lower signs refer to 
the /i and r flavors, respectively. The resonance terms, 
proportional to P^, appear where they should: in the 
neutrino sector for the NH, and in the antincutrino sec- 
tor for the IH. Note also that in the nonadiabatic limit, 
i.e., with set equal to one, the 2nd resonance be- 
comes irrelevant and the NH and IH cases properly re- 
duce to the same result. Note also that the asymmetry 
between neutrinos and antineutrinos is unsurprising: the 
solar medium contains no antimatter and in this sense 
provides a CPT- violating background [3lj |. 

In what follows, we will omit the small, explicit #13 
dependence from Eqs. (|24|) and (|25|) for brevity, but we 
will of course retain the primary #13 dependence in P c ' 1 
as given by Eq. (fT5|) . With #13 removed, the mixing is 
strictly tribimaximal, with its inherent ■<-> v T inter- 
change symmetry: the ± signs in Eqs. ([24]) and (1251) . 
which differentiate and v T , disappear. 

We now turn to the event types in the detector. A mi- 
nority of the events will be neutral current (NC) events. 
NC events yield no flavor information, and we ignore 
them. The majority of events will be charged-current 
(CC) events producing a charged lepton. Among these 
CC events, we note that the threshold in energy for pro- 
duction of a charged r is m r (1 + m T /2 tojv) = 3.47 GeV. 
Since this value falls in the middle of the sweet region, it 
presents a kinematic complication. It is therefore advan- 
tageous to consider observables that do not depend on 
CC r production. For his reason, we examine the ratios 

(W^/We), (W^/We), (We/We), and (Wp/Wp). The 

threshold for the muon CC is E — to m (1 + to m /2tojv) = 
110 MeV, well below the lower end (300 MeV) of our 
sweet region. For experiments with charge identifica- 
tion, e.g., the proposed large magnetized iron INO exper- 
iment, these ratios may be inferred from measurements 
of CC events, which necessarily contain a muon-track 
signature, as well as from v R CC events, which necessar- 
ily contain an electromagnetic plus hadronic shower and 
no muon track. For experiments without charge identi- 
fication, which can contain a much larger target mass, 
we must sum v and v and work with the single ratio 

U= (W^+W^/iWe+We). 

In practice, there will be experimental efficiencies that 
differ significantly for the detection of neutrinos versus 
antineutrinos, and muon versus electron CC reactions. 
Among other sources of these efficiency differences are 
the unequal CC cross sections for scattering of ver- 




FIG. 2: The normal-hierarchy ratio TZ NH = (W™ + 
W^ Ii )/(W^ lli + W^ H ) (top) and inverted-hierarchy ratio 

TZ m = {W™ + W™)/(W™ + W™) (bottom) versus the neu- 
trino energy E for the indicated #13 mixing angles. Here, 
Mwimp = 100 GeV, and the annihilation mode is assumed 
to be predominantly t + t~. The approximate location of 
the "sweet region," where our analytical results apply, is 
shaded. The star (★) indicates the location of the energy 
£1/2 at which the crossing probability is one half. The neu- 
trino mass-squared differences are 8m%\ — 8.0 x 10~ 5 eV 2 and 
dmh = ±3.0 x 10" 3 eV 2 , 02i and 832 are given their tribimax- 
imal values, and the CP-violating phase 5 is set to zero. The 
family of curves are quite sensitive to the value of 613, even 
below a degree. Relative changes in 1Z are 20% for the normal 
hierarchy (top) and 30% for the inverted hierarchy (bottom), 
with the centered energy E± / 2 showing quadratic sensitivity 
to 613 in agreement with Eq. (|20|) . We remark that the non- 
trivial behavior just below the sweet region results primarily 
from the lower-energy resonance as well as from the energy 
dependence of U m (E). 



sus v e versus versus v e . Over the energy range of the 
sweet region, the nature of neutrino scattering is quite 
energy dependent. Relevant cross sections include in- 
verse /3 decay, neutrino-electron scattering, quasi-elastic 
scattering, pion production, and deep inelastic scatter- 
ing. We note that efficiency differences may be advanta- 
geous if they can mitigate the mixing of v^, v e , D^, and 
v e signals. However, in this work we will simply assume 
uniform efficiencies for detection of v^, v e , D^, and v ei 
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FIG. 3: The normal-hierarchy ratio W^ H /W e NH versus the 
neutrino energy E for various #13 mixing angles. Input param- 
eters are identical to those in Fig. [5] The relative change in 
flavor ratio is about 40% across the transition region. Again, 
E1/2 (denoted by *) depends quadratically on #13, allowing 
a potential inference of this mixing angle. In the normal hi- 
erarchy, the 2nd resonance occurs in the neutrino sector and 
not in the antineutrino sector. 
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4: The inverted-hierarchy ratio W™ /W™ versus the 
neutrino energy E for various #13 mixing angles. The in- 
put parameters are identical to those in Fig. [2] The relative 
change in flavor ratio is almost a factor of two across the tran- 
sition region. Again, E^/i (denoted by *) depends quadrat- 
ically on 813 allowing a potential deduction of this mixing 
angle. In the inverted hierarchy, the 2nd resonance occurs in 
the antineutrino sector and not in the neutrino sector. 



rather than introduce another layer of complexity. 

The ratio K = (W M + W^)/(W e + W e ) is denned to 
weight equally the event numbers of neutrino and an- 
tineutrinos. Since the 2nd resonance occurs only in the 
neutrino sector (NH) or only in the antineutrino sec- 
tor (IH), but not in both, the contribution of the non- 
resonant sector (neutrino or antineutrino) will mitigate 
the contribution from the resonant sector. In this sense, 
1Z as defined is a conservative variable. 

The energy dependence of the ratio 1Z is presented 
without approximation in Fig. [3] One sees that even for 



this conservative variable, 20-30% changes in the relative 
flavor ratio occur across the adiabatic-to-nonadiabatic 
transition in the 0.3-10 GeV sweet region. Curves are 
parameterized by values of #13, with significant effects 
apparent all the way down to a fraction of a degree. 
The quadratic dependence of E1/2 (defined as the en- 
ergy where the nonadiabatic crossing probability is 50% 
and denoted in the figures by a *) on 9 13 is clearly seen. 
In favorable circumstances, this ultra-sensitivity of the 
spectrum to 613 could potentially be used to measure #13. 

Notable is that the change in the flavor ratio is about 
twice as large for the inverted neutrino mass hierarchy 
as it is for the normal hierarchy. We can dissect some 
of this behavior by turning to ratios defined only in the 
neutrino sector or only in the antineutrino sector. By 
excluding the non-resonant sector from the flavor ratio, 
we expect the change in the ratio across the adiabatic- 
to-nonadiabatic transition to be about twice that of the 
sector-summed result. That is, we expect ~ 40% ef- 
fects in the neutrino sector if the hierarchy is normal, 
and ~ 80% effects in the antineutrino sector if the hi- 
erarchy is inverted. In Figs. [3] and HI we see that this 
expectation is met. Consequently, if some experimental 
handle is available to even partially separate the neutrino 
and antineutrino data samples, then considerably greater 
discriminatory power becomes available. Put another 
way, across the sweet region, with the normal hierarchy, 
the neutrino flavor ratio should exhibit considerable en- 
ergy dependence and the antineutrino flavor ratio should 
not. Conversely, with the inverted hierarchy, the antineu- 
trino flavor ratio should show significant energy depen- 
dence and the neutrino flavor ratio should not. Thus, if 
neutrino-antineutrino discrimination becomes possible, 
then the observed energy dependence of neutrino versus 
antineutrino flavor ratios differentiates between the two 
possible mass hierarchies. 



B. Approximate results 

With the exact evolved flavor relations of Eqs. (|24p and 
in hand, we may take ratios and expand in powers 
of the small /f's given in Table |U To order f F , one finds 
for the normal hierarchy: 



n 



NH 




= 1 
= 1 

= 1 



27 
32 



£(i-j*)/ F + ^(i-j*)»/*, 



9 f 4- 27 f 
4 JF + — f- 



2 

F • 



F ■ 



(26) 
(27) 

(28) 

(29) 

(30) 
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These equations reflect the fact that with the nor- 
mal hierarchy, the 2nd resonance and the adiabatic-to- 
nonadiabatic transition lie in the neutrino sector and not 
in the antineutrino sector. 

The analogous calculation for the inverted hierarchy 
yields 



p=i 




tIH 



IH = x _ 


9 


2P c h ) f F - 




2P } C 1 ) 2 f 2 F 


,(31) 




0-/f 


f Q • fp , 






(32) 


D- 


9 


2P c h ) fp + 




2P h c f f F 


,(33) 


;).,. 


5(1- 

2 V 


2P c h )f F ^ 


■o-/l, 




(34) 


I).,- 


4 V 


IPc) fp + 


o • f F . 




(35) 



The equations here reflect the fact that with the in- 
verted hierarchy, the 2nd resonance and the adiabatic- 
to-nonadiabatic transition lie in the antineutrino sector 
and not in the neutrino sector. 

For comparison, we remind the reader that in the ab- 
sence of matter, the analogous vacuum-evolved ratios are 



(WJW e )„ m = (W /W r e )vacuum = 1 + f fp + I 



as 



has been established in Sec. IIVI For the value f F — 
0.18 appropriate to the r + r _ annihilation mode in the 
sweet region, this ratio is equal to 1.30. For the value 
fp = —0.09 appropriate to the bb annihilation mode in 
the sweet region, this ratio is equal to 0.9. Inserting ei- 
ther of these two fp values into the ratios calculated with 
solar-matter effects, Eqs. (|2"6"1) - ([33)) . gives very different 
results. 

We note that the change in flavor ratios over the tran- 
sition from adiabatic (P c ' 1 = 0) to nonadiabatic (P c ' 1 — > 
cos 2 #13 ~ 1) is twice as large for the antineutrino sector 
in the inverted hierarchy, as compared to the neutrino 
sector in the normal hierarchy. This factor of two is sub- 
tle. We reveal its origin in the Appendix. 

We may now use Eqs. ([2"6" |) - (|3"5)l to obtain the mag- 
nitude of changes in the flavor ratios as the neutrino 
energy E varies across the sweet region. The energy- 
dependent Pc(E) is a monotonic function and obeys 
P£(E -> 0) = and P' C }(E ->■ oo) = cos 2 <9 13 . Thus, 
at leading order in #13 the range of the crossing proba- 
bility is ^ Pc(E) ^ 1- We here list a sample of ratios 
of the quantities in Eqs. (|2"B"1) - (|3"5|) above and below the 
adiabatic-to-nonadiabatic transition, calculated to order 
fl, and then evaluated with the choices fp = +0.18 for 
the t + t~ mode and fp = —0.09 for the bb mode. In the 
normal hierarchy we have 



TpNH 

/c p=i 

nNH 



1 + ^/f + ^/I 



27 
32 



1 + 0.23 (r+r- 
1-0.09 (bb) 



, (36) 



(WNH/W e NH) P=0 



-1 + 9 f + 27 f 2 - J 1 + - 46 
- 1 + l fF+ 16 fF -\ 1-0.19 



while in the inverted hierarchy we have 



(37) 



n 



IH 

_P=1 



IH 

P=0 



1 + 7 h 



189 



0.50 (t+t- 
0.18 (bb) 



(38) 



(W™/W™) P=1 

(W™/Wf) P=0 



1 + 1.14 
1 - 0.32 



. (39) 

The value of these ratios of ratios in vacuum is unity, of 
course. The large flavor-ratio differences here are due to 
the adiabatic-to-nonadiabatic transition at the matter- 
induced 2nd resonance. For the normal hierarchy, the 
ratio of ratios for the neutrino changes by 23% and 46% in 
the t + t~ annihilation mode, and by 10% and 20% in the 
b b mode. For the inverted hierarchy, the ratio of ratios 
for the antineutrino changes by slightly more than twice 
that, by 50% and 114% in the t + t~ annihilation mode, 
and by 20% and 32% in the b b mode. The changes for the 
two hierarchies are of opposite sign. Together with the 
very little change expected in the W + W~ annihilation 
mode, one sees that just an inference of the existence 
and sign of the change can discriminate among the three 
most popular annihilation modes of solar WIMPs. 

The qualitative factors for the change in the ratio of 
ratios, given in Eqs. (|3"6"j) - (|3l?)) . are borne out in the more 
exact results shown in Figs. [2HU 



VI. DISCUSSION AND CONCLUSION 

We have shown that 

• For neutrinos from the t + t~ annihilation mode 
of solar WIMPs, 10-100% changes in neutrino flavor 
ratios at Earth are expected in the 0.3-10 GeV energy 
interval. This energy interval is a mapping across the 
adiabatic-to-nonadiabatic transition of the higher-energy 
2nd resonance in the Sun. Whether the change is 10% 
or 100%, or in between, depends on which flavor-ratio 
observable is available to the experimenter. 

• For the bb annihilation mode of the WIMPs, the flavor 
change is about half as large, and of opposite sign. 

• For the W + W~ annihilation mode, the flavor change 
is small. 

• Arising from the higher-energy resonance in the 
Sun, the effect occurs in the neutrino sector but not 
antineutrino sector if the neutrino mass hierarchy is 
normal, and in the antineutrino but not neutrino sector 
if the mass hierarchy is inverted. 

• Thus, a measurement of the sign and magnitude of 
this flavor change could indicate the WIMP annihilation 
mode, and a determination of the resonant sector, neu- 
trino or antineutrino, would indicate the mass hierarchy. 
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A plot of flavor versus energy determines the energy 
E1/2 at which the adiabatic-to-nonadiabatic transition 
probability is 50% and determines the width (in energy) 
AE of the transition. It turns out that 

• El/% s» AE, and each is proportional to (9f 3 . Thus, the 
system is over-constrained, which provides a clean signal 
for the underlying physics. 

• Moreover, the quadratic dependence on #13 provides a 
very sensitive probe for measuring the value of #13. 

• The adiabatic-to-nonadiabatic transition profile is 
measurable for # 13 as small as 0.5°. 

To establish experimentally a 30% change in the fla- 
vor ratios over the transition region with an n-sigma sig- 
nificance, the number of required CC v e and events 
is ~ 2 (n/30%) 2 . (The pre-factor 2 arises from the 
need for statistical accuracy on each side of the tran- 
sition.) Thus, for 3cr significance one needs roughly 200 
solar WIMP events separated from the atmospheric back- 
ground and spanning the sweet region between 0.3 and 
lOGeV. According to Eq. ©, a 200-event sample can 
be accumulated in the sweet region with a megaton de- 
tector in 1.5 years, if the WIMP annihilation cross sec- 
tion saturates its present experimental upper limit of 
3 x 10 _38 cm 2 . The annihilation cross section could have 
this value, or it could be orders of magnitude less. 

We end with mention of the many mega-detectors for 
low-energy neutrinos that have been proposed. They 
represent a presence of hope and vision that will shape 
the future. Proposed detector materials include water 
for Cerenkov signals (MEMPHYS, UNO, and HyperK), 
scintillators such as liquid argon (LENA, GLACIER, 
LArTPC, and TASD), and magnetic iron calorimeters 
(MIND, MONOLITH, and INO). The magnetized de- 
tectors can distinguish neutrino and antineutrino events, 
which enhances the signal explored herein by a factor of 
two, and thereby reduces the data required by a factor of 
four. But magnetized detectors cannot be built as large 
as unmagnetized detectors. 
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Appendix: The factor of two between antineutrino 
and neutrino adiabatic-to-nonadiabatic evolution 



To explain the relative factor of two at order fp be- 
tween the neutrino's and the antineutrino's adiabatic-to- 
nonadiabatic evolution, we use the density-matrix for- 
malism. In the flavor basis, the density matrix at Sun's 
core is p F = \ t + fp [2 \v T ){u T \ - |z^)(*Vl - \v e )(v e \]. 
Invoking the completeness relation ^ a |^ a )(^ a | = 1 
as well as maximal neutrino mixing to symmetrize 
(21^)^1-1^)^1) ->■ \{\y T ){v T \ + \v tl ){v ll \) shows 
that this matrix may be rewritten in terms of just the v e 
projector as pp -> | [ (1 + | f F ) 1 - § fp |^ e )(^e| ] ■ 

Now, with the normal hierarchy, and in the lower- 
energy adiabatic region, the v e emerges from the Sun as 
V3, and the and v T are left to emerge as v\ and 1/2 ■ In 
the higher-energy nonadiabatic region, the v e emerges as 
1/2 , while the and v T emerge as v\ and V3 . (Recall that 
we have chosen the sweet region to be below the nonadi- 
abatic onset E l NA of the lower-energy 1st resonance, and 
so the conversion across the 1st resonance is v e f2-) 

On the other hand, with the inverted hierarchy, and in 
the lower-energy adiabatic region, the v e again emerges 
as D3 and the P„ and v T emerge as and v^. However, 
in the higher-energy nonadiabatic region the v e emerges 
as v\ , while the and v T arc left to emerge as T)\ and vi . 

Thus, at lower energies we have v e — > v% and v e — > D3 
in the normal and inverted hierarchies, respectively, while 
at higher energies we have v e — ¥ v-i and v e — > v\ in 
the normal and inverted hierarchies, respectively. Re- 
placing the v e projector \v e )(v e \ in pp with the appro- 
priately evolved z/-mass projector or 1^2X^2 1 or 
I V3) (1S3 1 , and then calculating the final flavor- matrix ele- 
ments Wn = {vh\pf\ v ii) and W e = (v e \p F \v e ), one readily 
obtains the change in flavor ratios {W^/We) across the 
adiabatic-to-nonadiabatic transition. To order fp, the 
result is I fp (U e 2— U^+U^ — U^) for the normal hier- 
archy, and I fp (U e i — f/jui +Ufj,3 — U e s) for the inverted 
hierarchy. Putting into these expressions the values of 
tribimaximal mixing, one finds that the change is | fp 
for the neutrino sector in the normal hierarchy, and twice 
that for the antineutrino sector in the inverted hierarchy. 
Notice that this result requires the mixing values of v e 
and with all three neutrino mass states. In the end, 
it is a "conspiracy" of the tribimaximal values that gives 
the seemingly simple factor of 2. 
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